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effective Lagrangian approach. Cross sections and single-polarization observ- 
ables are investigated for the four charge channels, "fp — > ir + n, jn — > vr^p, 
7P — > 7r°p, and 7n — ► 7r°n. Leaving the ttNA coupling strength to be a free pa- 



rameter, we obtain a reasonably consistent description of these four channels 
■ from threshold to the first resonance region. Within this effective Lagrangian 

^sO ■ approach, strongly constrainted by the quark model, we consider the issue 

, of double-counting which may occur if additional i-channel contributions are 

O ! included. 
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I. INTRODUCTION 



Pion photoproduction has provided a wealth of information about baryon resonances. 
During the past three or four decades, extensive investigations have been carried out in both 
experiment and theory. In particular, the recent availability of high intensity electron and 
photon beams at JLab, ELSA, MAMI and ESRF has significantly improved the precision of 
pion photo- and electroproduction experiments. A large experimental database now exists, 
and a significant increase is expected once the current set of experiments has been analyzed. 

Pion photoproduction has been an important source, supplementary to itN scattering 
experiments, for establishing most of the well-known baryon resonances, while providing 
information on their photo-decay amplitudes. In the search for "missing resonances" , other 
meson production channels, to which these resonances might have stronger couplings, are 
now being extensively studied (See e.g. Ref. [I] and references therein). 

Apart from a few dominant states, a considerable model-dependence exists in resonance 
parameters extracted using phenomenological approaches to the data. This has complicated 
the comparison with resonance parameters derived from quark models. Historically, most 
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approaches have adopted a factorization of the meson interaction vertices, where the dy- 
namical information is absorbed into the resonance partial-decay widths and empirical form 
factors. Consequently, parameters for the meson-nucleon-resonance couplings and form fac- 
tors have been introduced. 

Such empirical schemes have been very important in analyses of data and the extrac- 
tion of resonance signals in pion photoproduction ||. A number of multipole fits, taking 
into account different dynamical aspects, are now underway. For instance, the unitary iso- 
bar model (MAID) ||, containing Born terms, 5 resonances and vector meson exchanges, 
succeeds in the description of data up to 1 GeV. Approaches adopting constraints from 
fixed-t dispersion relations are being re- visited and applied to the Delta region . Other 
approaches, using effective Lagrangians for the Delta resonance excitation and t-channel 
vector- meson exchange, can also be found in the literature [pHTTfl . The SAID fits fl2| . |T3l . 
based on a parametrization of different partial wave contributions, extend the analysis up to 
2 GeV. With some common features but quite different model constraints, these multipole 
fits hope to converge to a common result and obtain, as near as possible, model-independent 
information on the resonance excitations. 

There is a clear need to treat all resonances consistently, and to understand the rela- 
tion between the s- and ^-channel resonances and t-channel meson exchanges. A recently 
developed quark model framework |Q, augmented by an effective Lagrangian approach to 
reaction dynamics, provides a good starting point. The main feature of this model is the 
introduction of an effective chiral Lagrangian for the quark-pseudoscalar-meson coupling in 
a constituent quark model. Unlike most previous quark models, which were generally based 
on factorization of the strong interaction vertices, the pion is treated as an elementary par- 
ticle. As a result, one can explicitly calculate the tree level diagrams for pion production 
reactions. Here, the quark model wavefunctions for the nucleons and baryon resonances 
provide a form factor for each interaction vertex, and all the s- and ^-channel resonances 
can be consistently included. 

This model has the advantage of being able to describe a large photoproduction database, 
employing only a very limited number of parameters within a microscopic framework. Appli- 
cations of this model to the r] |fT5-I8|| and K flT9| , |20|| meson photoproduction have been quite 
successful, and this has motivated our study of the very extensive pion photoproduction 
database. 

The quark model's well-known underestimation of the electromagnetic (EM) transition 
amplitude for the Delta resonance makes this resonance region particularly interesting. As 
suggested in Ref . |Tl| , the "bare" jN — > A vertex could be more directly related to the quan- 
tity given by the quark model derivation of the Delta EM transition. A direct examination 
of the Delta excitation in — ► A — ► nN might shed some light on this question. 

This paper presents both quantitative and qualitative investigations of pion photopro- 
duction. The challenge to describe the N* resonance excitation with explicit quark and 
gluon degrees of freedom is by no means trivial, since the correct off-shell behavior of those 
intermediate resonances is required. Also, a clear definition of the nucleon Born terms, 
associated with the gauge invariance requirement, is essential for this effective theory. In 
this study, we concentrate on the energy region corresponding to E 1 < 700 MeV, where the 
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role played by the Born terms [] and the low-lying resonances, in particular the A(1232), 
Sii(1535) and Di 3 (1520), can be clarified. Qualitative tests have been made in order to 
compare this model to a typical isobaric approach. Here we consider the role played by 
the t-channel vector meson exchanges in isobaric models, and the effect of neglecting the 
w-channel resonance contributions. 

In Section II, we outline the formalism aspects of our approach. In Section III, results for 
cross sections and single polarization asymmetries for the four charge channels, jp — > n+n, 
7P — > 7t°p, jn — » tt~p, and 772 — > 7r°n will be presented. The role of the t-channel vector 
meson exchange will also be discussed. Conclusions are drawn in Section IV. 



II. THE MODEL 

Before we begin a detailed analysis, a brief review of this model is necessary. 



A. The effective Lagrangian 



For pion photoproduction, the low energy theorem (LET) provides a crucial test near 
threshold. As shown in previous investigation by Li [T2TJ , to recover the LET, one has to 
rely on the low energy QCD Lagrangian which keeps the meson-baryon interaction invariant 
under the chiral transformation. Combining the low energy QCD Lagrangian with the quark 
model, we introduce the quark-meson interaction through the effective Lagrangian (TJ : 



C = ipMid" +V + 75^) - m]ip + 
where the vector and axial currents are 



(1) 



(2) 



and the chiral transformation is, 



(3) 



where f m is the decay constant of the meson. The quark field ip in the SU(3) symmetry is 



if, 



(4) 



and the meson field m is a 3®3 matrix: 



We use "Born terms" here to denote the amplitudes from a Born approximation, in which the 
nucleon pole terms, pion pole and contact term are included. In the following Sections, we use 
"nucleon pole terms" to denote the s- and w-channel nucleon exchange amplitudes. 
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(5) 



where the pseudoscalar mesons n, r\ and K are treated as Goldstone bosons. Thus, the 
Lagrangian in Eq. (HD is invariant under the chiral transformation. Expanding the nonlinear 
field £ in Eq. (|3]) in terms of the Goldstone boson field m , i.e. £ = 1 + i(p m / fm + • • •, we 
obtain the standard quark-meson pseudovector coupling at tree level: 



j Jm 



(6) 



where ipj (ipj) represents the jth quark (anti-quark) field in the nucleon. 

It is still not clear whether the Goldstone bosons couple to the nucleon through a pseu- 
doscalar or pseudovector coupling, or even both. To our present knowledge, at low energies, 
the pseudovector coupling satisfies partial conservation of axial current (PCAC) and is con- 
sistent with the LET and chiral perturbation theory to leading order, while the high energy 
study prefers a pseudoscalar coupling. As pointed out in Ref. |22[], the operators for the pseu- 
doscalar and pseudovector coupling have the same leading order expression at quark tree 
level. Therefore, Eq. (^) can be regarded as a reasonable starting point for investigations of 
pion photoproduction in the resonance region. 

The quark-photon electromagnetic coupling is 



^e = -E e ,7^(k,r), 



(?) 



where the photon has three-momentum k, and the constituent quark carries a charge Cj. 
The photoproduction amplitudes can be expressed in terms of the Mandelstam variables, 



M 



fi 



M% + M s fl + M« + M%, 



(8) 



where M^f is the seagull term and Mj iy and Mj i represent the s-, u-, and t-channel 
processes as illustrated in Fig. [I]. As shown in Ref. JTJJ, the seagull term is composed of two 
parts, 



M)\ = (NflH^Ni) +i{N f \[g e ,H m ]\N i ), 
where |A^) and \Nt) are the initial and final state nucleon, respectively, and 



j J in 



is the contact term from the pseudovector coupling, and 



ikr, 



(9) 



(10) 



(11) 



comes from the transformation of the electromagnetic interaction in the s- and u- 



channel [23.14 
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The s- and w-channel amplitudes have the following expression: 



where 



^2(Nf\H m \Nj){Nj 



E; 



OJ, 



-h e \Ni) 



+ iuj. 



E; 



LOr. 



E 3 

Nj){Nj\H m \Ni), 



h e = ejTj ■ e 7 (l — cx.j ■ k)e 



ikr, 



(12) 



(13) 



and k = k/u; 7 is the unit vector in the direction of the photon momentum. 



The nonrelativistic expansions of Eqs. ( |13l) and (y) become [14 



E 



6 1 rr 

2m,- 3 



(e 7 x k) 



and 



HI 



E 



-o\--P f + 



E f + M f 3 1 Ei + Mi 3 



a-j Pi - a-j ■ q + 



ikr, 



°i • Pi 



J 



e ~ iq ' rj 



(14) 



(15) 



where Mj (M/), £7j (-E 1 /) and Pj (P/) are mass, energy and three- vector momentum for the 
initial (final) nucleon; and pj are the internal coordinate and momentum for the j'th quark 
in the nucleon rest system. Note that, gA, the axial vector coupling, relates the hadronic 
operator cr to the quark operator <Tj for the jth quark, and is defined by, 

(NflY.I^m^gAiNflcrm. (16) 

j 

The transition amplitudes of pseudoscalar meson photoproduction can generally be ex- 
pressed in terms of standard CGLN amplitudes P3fl , i.e. 

Mfi = J • e 7 , (17) 

where J is the interaction current and can be related to the CGLN amplitudes /i,2,3,4 : 



(cr • q)(k x cr) cr ■ k 
fiv + 1/2- " + /s^^q 



q k 



|q| |k| 



, <* • q 
U — 5~q- 
q 2 



(18) 



Alternatively, one can express the transition amplitudes in the helicity space in terms of the 
T matrix: 



Hi = 




+1/2|T|A 7 


= +l,Ai = 


-1/2) 


H 2 = 


(A/ = 


+1/2|T|A 7 


= +l,Ai = 


+1/2) 


H 3 = 


(A/ = 


-1/2|T|A 7 


= +1,A; = 


-1/2) 


H A = 


(A/ = 


-1/2|T|A 7 


= +l,Ai = 


+1/2) 



(19) 

where \ and A f are helicities of the initial and final nucleons and A 7 is the helicity of the 
photon. Amplitudes with A 7 = — 1 are not independent of those with A 7 = +1 due to 
parity conservation. The CGLN and helicity amplitudes may be related through a standard 



transformation 25 



5 



B. Transition amplitudes in the harmonic oscillator basis 



The seagull term in this model differs from the traditional definition due to the appear- 
ance of a transformed electromagnetic interaction coupling to the meson at the same vertex. 
This term can be worked out explicitly in the SU(6)®0(3) symmetry limit: 



-(k-q) 2 /6a 2 



+ 



Ei + Mi Ef + M 



a ■ e 
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(20) 



where the exponential factor is the corresponding quark model form factor in the harmonic 
oscillator basis. 

The t-channel charged pion exchange amplitude can be derived by treating the exchanged 
pion as an elementary particle: 



M% = e 



-(k-q) 2 /6a 2 e m (Mf + Mi) 

q ■ k 



a ■ q 



a ■ k 



E f + M f Ei + Mi 



q- e 



7' 



(21) 



where q and k are four-vector-momenta of the pion and photon, respectively. 

As illustrated in Eqs. ( |2"0"| ) and (|2"T|), the leading order amplitudes from chiral perturbation 
theory are reproduced. The quark model modifications to these two terms come from three- 
body effects, which add an additional term (the second term in Eq. (PD|)) to the amplitudes. 
Note also the appearance of a form factor, which is essential to sustain the forward "spike" 
in 7r + production. 

Generalized expressions for the s- and n-channel amplitudes are: 



Ml = (M 2 * + M 3 S K 



-(k 2 +q 2 )/6a 2 



(22) 



and 



M 



fi 



(M% + M3 U )e~ (k2+q2)/6a2 



(23) 



where the M 3 and M 2 represent transitions in which the photon and meson couple to the same 



quark or different quarks, respectively. The general framework was presented in Ref. [14 



Here, we present the transition amplitudes in terms of the harmonic oscillator shell n as 
follows: 



Mi 1 

— = — [ig v K • (e 7 x k) - cr • (A s x (e 7 x k))] 
93 2m g 



M„ 



+ 



^7^m n _|_ u h 

H q 2m q 



)<r ■ e 7 + 



2u. 



+ T7T of* ■ ke 7 ■ q- 



M„ 



cr ■ A s e 7 ■ q 



n\(Pi ■ k - nMu h ) 

M n 

(n 



k ■ q N 

3^ 



n - 2)\(Pi ■ k - nMoo h ) \3a 2 



M 

l)!(Pi • k - nMuj h ) \3a\ 

n-2 



n—1 



(24) 



and 



6 



M|(-2) 
9 S 2 



2m, 



WvK ■ (e 7 x k) - g' a a ■ (A, x (e 7 x k))] 



M n /k_q N 

n\(Pi ■ k - nMuJh) \ 3a 2 



+ 



u; 7 u; m 



1 + 



/ W 7 \ 



cr ■ A s e 7 • q 



(n-l)!(P • k-nMu h ) \ 3a 2 



n-l 



+ ttt ^cr • ke 7 • q- 



k • q 



n-2 



18n q a 2 "~ 7 ^(n-2)!(Pi-A;-nM^) V3a 2 
result for the s-channel, while 



(25) 



M 3 n 1 

1 ' 

~ 6 



[ig v A u • (e 7 x k) + cr • (A u x (e 7 x k))] 

2u. 



M„ 



U-yU m Q _|_ W 7 



0g 

U^U m 

18[i q a 2 

and 
M 2 u (-2) n 



2m/ CT ' e - + 



-cr A„e 7 • q 



n\(P f ■ k + nMu h ) 
M n 



k • q N 

Ito 2 



k • q 



n— 1 



cr • ke 7 ■ q- 



M n 



92 



2m, 



k-q 

n - 2)\(P f ■ k + nMu h ) \3a 2 



[ig' v A u ■ (e 7 x k) - g' a cr ■ (A u x (e 7 x k))] 



n- l)\(P f ■ k + nMuOh) \3a 2 

n-2 



(26) 



Mr. 



n 



/k 1 q X 

(Pfk + nMu h ) \ 3a 2 t 



1 



^■yUJrn , UJ-y . 



6 [ //g 



'2m n 



M n 



2^ 

7+ a 2 

k • q 



cr A M e 7 • q 



n-l 



n 



1)!(P, -k + nMuJh) V 3a 2 



u; 7 u; m 
18/i g a 2 



cr • ke 7 ■ q- 



M n 



k • q 



n-2 



-2)!(P/ • k + nMu h ) \ 3 



ex 



(27) 



are corresponding terms for the it-channel. Vectors A s and A u are determined by the meson 
transitions in the s- and it-channels: 



A s 



E f + M 



+ 1 q, 



/ 



(28) 



and 



A.q, 



U r . 



+ 



UJ r . 



Ei + Mi E f + M 



E f + Mj 



+ 1 q 



(29) 



In Eqs. (|2^)-(|27j), Pj and P/ are four- vector momenta of the initial and final state nucleons in 
the total cm. system; M n is the mass of excited state in the nth shell, while Uh (= a 2 /m q ) 
is the typical energy of the harmonic oscillator potential. The factor M n j (Pj ■ k — nMujh) 
and M n j (P/ • k + nMuh) have clear physical meanings in the 5- and it-channels that can be 
related to the propagators. 
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The quark level operators have been related to the hadronic level ones through ^-factors 
defined as below: 



23 u H^/l£^Vil^>M, (so) 



92 = (N f \ y £e j i i at\N i )/g A , (31) 



9v = -^-(Nf\ y £e j i j m i (32) 
93 9a ,■ 



9v = ^r( N f\ E e i^ • <^>> (33) 



and 



9a = 77-(JV/EeA<T ( x ffi ),|JVi). (34) 
4*72 i^j 

Numerical values for these g-factors can be explicitly calculated in the SU(6)®0(3) symme- 
try limit M . 



So far, the resonance contributions have not been explicitly separated out. The in- 
termediate states are still degenerate in the quantum number of the harmonic oscillator 
shell. Notice that the factor M n /(Pj • k — nMujh) can be written as 2M n /(s — M 2 ), where 
s = W 2 = (Pi + k) 2 is the square of the total cm. energy, we thus substitute a Breit-Wigner 
distribution for the resonances, 2M R /(s - M\ + iM R T R (q)), where the resonance mass and 
width effects are taken into account. Explicitly, all the contributing resonances with n < 2 
in the quark model symmetry limit can be included. In pion production, this is the place 
where the imaginary part of the transition amplitude comes out. The role of the imaginary 
part can be investigated more directly in the polarization observables, e.g. polarized target 
asymmetry T and recoil polarization asymmetry P, where direct interferences between the 
real and imaginary part are highlighted. 

Consequently, we must separate out the resonance excitation amplitudes for each n. For 
n — 0, the contributing terms are the Delta resonance excitation and the nucleon pole terms. 
One can see that only the first terms in Eqs. (p4])-(p7|) can contribute. For the s-channel, 
we have 

M s (n = 0) = -JL[i(g> 9v + g y v )A s ■ (e 7 x k) 

-G?3 + 9l9> ■ (A s x (e 7 x m Pt , k M l M2/2 e~ (k2W)/ea \ (35) 

where SM 2 /2 denotes the mass square difference between the intermediate state and initial 
state nucleon. The amplitude for the spin 1/2 nucleon pole term is, 



S 



M s (nucleon) = (N f \H m \N(J = 1/2))(N(J = l/2)\h e \Ni) 



- -ItLv . A cr ■ (e x kl 2Mn r -(k 2 +g 2 )/6a 2 ( o 6) 
~ 2m* (C7><kj S -M2 e ' W 

where we have used 5M 2 /2 = and Pi ■ k = (s — M%)/2 for the nucleon pole; \ii is the 
magnetic moment of the initial nucleon in terms of the proton magnetic moment e/2m q . In 
this way, the Delta resonance excitation amplitude is derived, 

M S (A) = M s (n = 0) - M s (nucleon) 

= ~^~~[^3^ + 929v ~ Vi) A s ■ ( e 7 x k ) 

-(^ + 9t9' n - tk)<r ■ (A s x (e- x k))l ^ e -(k^W)/6a 2 ( 37) 

where M /(Pi ■ k - 5M 2 /2) = M A /(s - M 2 N - (Ml - M 2 N ))/2 is used and the Breit-Wigner 
distribution is introduced after the separation of the spin operators. 

Similarly, the Delta resonance and nucleon pole terms in the it-channel with n = can 
be separated: 

M u (nucleon) = (N f \h e \N(J = 1/2)}(N(J = l/2)\H m \Ni) 

^La ■ (e 7 x k)<r • A u -%e-( k2+ ^/ fe2 , (38) 

and 



2m q ^ ' U u-M 2 N 



M U (A) = -J-[i(g-g v + g u 2 g' v - ^ f )A u ■ (e y X k) 

+(9 U s - 9Wa - A*/)o- • (A u x (e 7 x m^^f^ W)/ " a \ (39) 

where [if is the magnetic moment of the final state nucleon in terms of the proton magnetic 
moment e/2m q . 

Several points can be learned from Eqs. (|2"4[)-(|27D. First, the nucleon pole terms and Delta 
resonance transition only involve the cm. part of their spatial wavefunctions. Therefore, 
only the first terms in these equations contribute to the amplitudes. For resonances with 
n > 0, the internal motion of constituents will be involved. Specifically, terms relating to 
(n — 1) are due to correlations between cm. motion and internal ones, while (n — 2) terms 
are due to correlations between internal motions at two vertices. 

Secondly, amplitudes for processes having the photon and meson coupled to different 
quarks are relatively suppressed. This can be seen clearly through the factors (— 2) n . In 
Ref. this qualitative feature was discussed, but not shown explicitly. Here, we show 
how the indirect diagram can be exactly calculated, and show that the direct diagram will 
become dominant with increasing energy and the excitation of higher states. 

Notice that in the degeneracy limit, the sum over n in Eqs. (p4|)-(p7j) gives: 



= (O s + O u )e-^ 2 / 6a2 , (40) 
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where O s and O u represent operators independent of n, and recovering the factor e~( k ~ q ) 2 / 6a2 
is essential for the theory to be gauge invariant. Although like many other phenomenological 
approaches this model does not have unitarity, such a form factor prevents certain terms from 
violating unitarity. One can see that at high energies the degeneracy limit can be recognized 
by the dominance of the direct diagram. How to restore the unitarity in a general framework 
should be a next step of this investigation. 

From Eq. (|37|), the analytical expression for the Delta multipole can be derived, 



/\/f 3 /2 1 



' jJi " + 1 



2m q [E f + M f 



2M A 3 -(k 2 +q 2 )/W 



Mi + iM A T A 



e-^ k+q ^ to , (41) 



where qr = g^gv + g^g'v ~ an d g-wNN has been taken into account. The real and imaginary 
parts of the Delta multipole iWf{. 2 are calculated and shown in Fig. |^. We shall discuss the 
quark model form factor effects in the following Section. Therefore, it would be useful to 
present the calculation of M 1+ without the exponent, which comes from the spatial integral 
and serves as a form factor for the multipole. As illustrated by the dotted and dot-dashed 
curve in Fig. |2], apparent deviations from the experimental data occur with the increasing 
energies. In another word, the quark model turns to be indispensible to account for the 
correct energy dependence. 

Multipole E 1+ vanishes in this approach for the s-channel Delta resonance. Experimen- 
tally, Ei+ is found to be much smaller than M\+ |13|J4]]. 



III. RESULTS AND ANALYSIS 

In this Section, we present our study of reactions, yp — > n + n, yp — > n p, •yn — > ir~p, and 
7n — > 7T°n with the same set of coherent parameters. The Goldberg- Treiman relation, 

9aM n 

g^NN = —f. — , (42) 

J TV 

relates the axial vector coupling g& to the well-known g n NN coupling, where /„- is the pion 
decay constant. Note that gA in this model is an overall constant, and can be calculated in 
the quark model. However, the quark model predicts a rather large value: \gA\ = 5/3 for 
charged pions and 5/3\/2 for neutral pions. Consequently, g^NN given by Eq. (^2|) is not 
a good input for our purpose. On the other hand, g^NN is a well-determined number, and 
we shall therefore "fix" the "parameter" g^NN = 13.48 in our calculations. This is a strong 
constraint on the theory, as the Born terms (including seagull term, t-channel pion exchange, 
and nucleon pole terms) are completely fixed. We shall follow the quark model predictions 
for relative signs and strengths in order to study the four charge channels coherently. 



A. The charged pion production reaction 

A distinguishing feature of tt + photoproduction is the forward "spike" and the dip near 
—t = mjj. seen in its differential cross sections. Multipole analyses based on a unitary isobaric 
model and studies utilizing fixed-t dispersion relations suggest that this feature is related 
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to an interference between the Born terms and the Delta resonance "excitation" (in the 
naive quark model, the Delta resonance is the ground state for isospin 3/2 baryons). The 

3 /2 

M 1+ multipole for the Delta resonance dominates the cross sections and single polarization 
asymmetries over a wide energy range. Reproducing this combination of features is non- 
trivial in a microscopic model. 

In Fig. H we show the calculations for reaction 7^ — > TT + n in the SU(6)®0(3) symmetry 
limit. The dotted curves are calculated in this limit, having only an overall quark- meson cou- 
pling parameter, which is in fact not "free" for the Born terms due to the Goldberg- Treiman 
relation. Near threshold, the cross section is reproduced reasonably, which is consistent 
with the leading order calculation of chiral perturbation theory. The dip at —t = m% clearly 
originates from the interference between the t-channel pion exchange and the seagull term. 
At E 1 < 220 MeV, the Delta resonance has only small interference. 

Interference from the nucleon pole terms becomes important from E 1 = 220 to 260 MeV. 
At the lower energy limit, the i-channel pion exchange and the seagull term dominate over 
other processes, while for higher energies, the Delta resonance dominates. The influence of 
the nucleon pole terms can be seen clearly in the polarized beam asymmetries. As shown 
by the dotted curve in Fig. |4] at E 1 = 260 MeV, for which the nucleon pole terms are elim- 
inated, the interference from the nucleon pole terms generally produces large asymmetries 
at intermediate angles. 

With increasing energy, we find that the cross sections are underestimated over the 
Delta resonance region in the SU(6)®0(3) symmetry limit. This suggests a failure of the 
symmetric quark model for the Delta resonance. However, notice that the "dip" feature 
is still sustained over this region, and we assume that the Delta excitation has a "good" 
form factor from the quark model. Thus, we empirically treat the 7riVA coupling strength 
as a free parameter, which will be fixed by the experimental data. The solid curves shown 
in Fig. |3] denote the calculations with Ca = 1.7474, where Ca = 1 is the strength in 
the quark model symmetry limit. The enhancement of the Delta contribution significantly 
improves the description of the experimental data. Compared with the dotted curves, the 
cross sections at the extremely-forward and intermediate angles are both enhanced. 

The differential cross section for the 7r + production changes rapidly from threshold to 
the Delta resonance region. After that, it remains stable up to E 7 ~ 700 MeV, where 
resonances of the second resonance region start to interfere. The challenge for a microscopic 
approach is not only to reproduce the dramatic changes at low energies, but also to sustain 
the forward peaking to high energies. It is quite natural for us to achieve the first goal 
in this model. That is, the enhanced Delta resonance succeeds in reproducing the drastic 
change of the cross sections at the first resonance region. For the second goal, we find that 
with only one parameter, the strong forward peaking can only be sustained up to E 1 « 500 
MeV. This result has non-trivial implications. It suggests that the quark model within an 
effective Lagrangian provides correct signs and even reasonable form factors for the Delta 
excitation and nucleon pole terms. As illustrated in Section II, the Delta excitation and the 
nucleon pole terms have simple structures coming from the harmonic oscillator shell n = 0. 
Clearly, a self-consistent treatment [36j of these ingredients is essential to any viable model. 
We suggest that the tree level calculation, based on the quark model, may have included 
the main ingredients, (e.g. relative signs and form factors), even though its description of 
the nucleon pole and Delta resonance is very simple. 
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An interesting question arising in this work is the role played by the ^-channel resonance 
contributions. Generally, this part has been neglected in isobaric models, nor is it included 



in traditional quark model calculations, due to empirical considerations In the present 
calculation, we find that the u-channel process tends to decrease the forward peaking. In 
Fig. |3], the dashed curves denote calculations with the tt-channel of n > 1 neglected, which 
enhances the forward peak above the Delta resonance region. Since the full calculation 
underestimates the forward peaking slightly, the neglect of the u-channel of n > 1 seems to 
follow the data more closely. This feature seems consistent with findings of Ref. ||. There, 
the tt-channel resonance contributions were neglected, and an overall strong forward peaking 
was observed. 

Polarization observables are sensitive to resonance contributions, providing a possible 
way to clarify the role played by the Delta resonance. In Fig. |], the polarized beam asym- 
metry E is calculated for 8 energy bins. The results are generally in agreement with the data 
at E 1 < 300 MeV. However, some discrepancies are found at = 350 and 400 MeV, which 
are sensitive to the w-channel nucleon pole rather than the w-channel resonances (n > 1). As 
shown by the dashed curves, neglecting the n > 1 ^-channel resonance does not change the 
solid curves significantly. At « 700 MeV, the Sii(1535) becomes a strongly interfering 
source. The enhancement of this asymmetry at 9 = 130-140° is evidence for the existence 
of the Sii(1535) resonance. 

The presence of the 5 , n(1535) as a state of representation [70, 2 8] in the quark model 
accounts for E naturally up to 750 MeV. Compared with the precise measurement of 
GRAAL (44| , we cannot produce the structure observed at backward angles above 800 MeV. 
As suggested by the isobaric approach ||, a small Sii(1650) contribution can reproduce the 
data reasonably. In our model calculation, the Sii(1650) is absent in the proton reaction due 
to the Moorhouse selection rule ]37j . The breaking of the symmetric quark model will intro- 
duce mixing between states of different representations, e.g. the Sii(1535) and 5ii(1650). A 
more realistic model taking into account such a mixing mechanism is clearly required above 
the second resonance region. The Sii(1650) has large branching fraction to ttN states ]38] 



Calculations for the polarized target asymmetries are presented in Fig. g, and compared 
with existing experimental measurements. Given that only one parameter has been intro- 
duced, the results should be regarded as consistent with the data from threshold to E 1 « 
500 MeV. At E 1 = 220 MeV, our results underestimate the data, however this feature is 
consistent with the SAID fit [fTSf . 

Calculations for the recoil polarization asymmetries are presented in Fig. ^, which are 
consistent with the data in the first resonance region. 

The reaction — > ir~p is calculated using the same set of parameters determined 
in the 7r + production. We present the results for the differential cross sections in Fig. [7]. 
Although large uncertainties exist within the data, our calculation is in good agreement 
with experiment up to E 1 « 400 MeV. Interestingly, these results, which can be regarded 
as predictions of this approach, are very close to the analyses of Ref. Q . It is worth noting 
that similar structures as found in jp — ► 7i + n (the "dip" and "spike") are also present here, 
and due to the same mechanism. 
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B. The neutral pion production reaction 



In comparison with the charged pion production, the neutral-pion channels are relatively 
simple in this model. The contact term and the ^-channel pion exchange are eliminated in the 
effective interaction since these amplitudes are proportional to the charge of the produced 
meson. The nucleon pole terms and the Delta excitation therefore dominate over other 
processes near threshold. 

In Fig. [8], the differential cross sections for 7^ — > 7c°p are presented at several energies. 
In the SU(6)cg>0(3) symmetry limit, as shown by the dashed curves, the cross sections are 
underestimated by the Delta excitation. Similar to the feature arising from the charged 
pion production channels, we need to enhance the Delta excitation strength to reproduce 
the data. 

The calculations of the single polarization asymmetries are presented in Fig. [11] and 
compared with data. 

In Fig. [12|, the calculated cross sections for 772 — > TT°n are presented. So far, there are 
only sparse data available for this channel. 

The above results for the cross section and single polarization asymmetries suggest an 
overall quantitative agreement with the data from threshold to the first resonance region, 
while qualitatively, the data up to the second resonance region can be explained. The lack 
of quantitative agreement above the Delta resonance region was expected, given that only a 
minimum number of free parameters are used here. However, by using only a minimal model, 
it has been easier to identify key ingredients responsible for those trends we do reproduce. 



C. Quark model form factor and the helicity basis 

As mentioned previously, in 7r + photoproduction, the most prominent features seen in 
the cross section are forward peaking and the dip at —t = m^, which is attributed to the 
Born terms. Our results also reproduce this feature. Some new ingredients appearing in 
this approach concern the roles played by the Born terms and the Delta resonance, and the 
influence of their associated form factors. 

As found in previous studies, the Born terms deviate significantly from the experimental 
data at intermediate and backward angles as photon energies increase to the GeV level. 
The cross section due to Born terms alone is much larger than the data suggest. One 
possible explanation is that the Born terms are cancelled by resonance contributions away 
from the forward peak. As discussed by Barbour, Malone and Moorhouse in a fixed-/; 
dispersion relation |p6|| , the real parts of the resonance amplitudes tend to cancel the Born 



terms at —t > m^, while the region —t < is slightly enhanced by low- lying resonance 
contributions. 

In Fig. [TBI, we illustrate the results for the Born terms and Born terms plus Delta exci- 



tation, with and without the quark model form factors, respectively. Clearly, form factors 
are vital in the quark-model description, though no free parameters have been introduced. 
Comparing the full result to one in which the form factors are switched off, we see potential 
problems for those who compare quark-model results directly to fits (such as SAID and 
MAID) which do not introduce form factors. An interesting extension of this work would be 
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the consideration of forward cross sections at higher energies, where the influence of form 
factors is nebulous [j8l| , p6 . 

To end this Section, we present a comparison of energy evolution of the Born terms plus 
Delta helicity amplitudes calculated by this model with a SAID analysis [jS2" |. The four 
independent helicity amplitudes are calculated following the convention of Ref. M, 



da 



l|q 
2 Ik 



4 

E 

i 



(43) 



where 8 cm is the angle between the incoming photon momentum k and outgoing meson 
momentum q in the cm. system. In Fig. |14|, the helicity amplitudes are presented at five 
angles, 9 cm = 0°, 45°, 90°, 135°, and 180°. At 9 cm = 0°, only H 2 has non-zero values, while 
all the other elements vanish. In the backward direction, the no n- vanishing element is H^. 
Compared with the SAID analysis, an overall agreement is obtained up to 500 MeV. 



D. i-channel vector meson exchange 

A long-standing question concerns the role played by vector-meson exchange in low- 
energy pion photoproduction. According to the duality argument |83|, the introduction of 



vector meson exchanges, along with a complete set of s-channel resonances, might result 
in double counting. In practice, a systematic inclusion of all s-channel resonances at the 
hadronic level is not available. Empirically, t-channel vector-meson exchange may account for 
incomplete s-channel resonance contributions, which however makes the duality hypothesis 
more ambiguous. 

Given the results presented in the previous Subsections, the quark model framework, 
with an effective Lagrangian, could address this question in pion photoproduction from a 
more fundamental level. As seen in the cross sections for tt + production up to 700 MeV, 
forward peaking above the Delta resonance has been successfully sustained up to 500 MeV. 
This could reasonably illustrate that the effective Lagrangian has been sufficient to describe 
the data over the first resonance region. In order to consider the effect of possible double 
counting between the i-channel vector meson exchange and s- and it-channel resonances, 
we compare models including various subsets of these diagrams. Our purpose is to clarify 
whether the behavior of those higher excited states (terms) would be similar to the inclusion 
of vector meson exchange, particularly to compare with an 'isobaric' model, where the u- 
channel resonance contributions are neglected. 

We shall introduce the following effective Lagrangians for vector-meson exchange: 



£ y = e i0^L £ d c A P d i V 5 7r, (44) 



for •jtiV coupling and 



C Vqq = 9v qq ^ 1, + 7^-^ V^, (45) 



\ 2m q 



for the quark-vector-meson (V-qq) coupling; A and V denote the photon and vector meson; 
ii denotes the pion field and -0 (-0) denotes the quark (antiquark) field; g 17T y an d gv qq are 
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coupling constants. Note that we treat the V-qq coupling at quark level in order to be 
consistent with our framework. In this way, there is once again no need to introduce free 
parameters for the vertex form factors. In addition, a simple current analogy will relate the 
Qvqq to 9vnn, constraining this term as well. 

Some simple algebra gives the transition amplitude, 



9V-K-y9Vqq e y ' ' 

m n (t - mf) 
Kn ( ml 



9t 



1 + 



Ef + M f 
1 



+ 



Bk + Mi 



+ 



2m q \Ef + M f \Ef + M f E{ + M, 



-k ■ q 



q 



fk x e„ 



+9a 



w 7 q 
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u; m k 2 



E f + M f Ei + M, 



2m„ 



u m k + cu 7 q 



E f + Mf 
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E f + M 
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ier ■ kq • e n 



(46) 



where k ■ q 



k • q is the four-momentum product; The exponent comes from the 



nucleon wave functions, which plays the role of a form factor; is the axial vector coupling 
and defined in the quark model as Eq. fllBD , i.e. 



{Nf\Y J i]<T i m=9A{Nf\<T\N i ), 



(47) 



where JJ is the isospin operator for the exchanged vector meson. The other factor g t comes 
from the isospin space, 



g t = (N f \Y,im). 



(48) 



Analogy between the quark level operator and V-NN coupling gives: 



9t9Vqq — 9VNN 

9v qq/ , . \ 9vnn , 

9A (1 + Kq) = I 

m q m N 



1 + K V ) 



(49) 



where m q = 330 MeV is the constituent quark mass. In Ref. [84J], a similar relation was 
investigated, but only the vector current was introduced for the V-qq coupling. We shall 
use the commonly used values for gvNN and Ky to constrain the values for g Vqq and K q . 

In the 7r + production, we adopt the values g P NN = 3 and n p = 3.71 as inputs. With the 
quark model values g& = 5/3 and g t = 1 for 7p — > n + n, we derive g pqq = 3 and n p = —0.0064. 



In 7p — > tc°p, we adopt g t 



ujNN 



9, and k u 



-0.12. With g% = 1 and gf = 3, g mq = 3 and 
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k% = 0.2 are derived for the u exchange, and with g p A = 5/3\r2 and gl = l/y/2, g pqq = 3 
and k 9 = —1.99 are derived. 

In Figs. ^ and IB, we show the calculations of observables (da/dQ, E, T and P) with the 
t-channel vector meson exchange (VME) for 'jp — > 7c + n and 7p — > ir p, respectively. Three 
energy bins are investigated. We use s+u to denote the effective Lagrangian calculations, 
while s+t denotes calculations suppressing the u-channel resonance but including the t- 
channel VME. Finally, we use s+u+t to represent the full calculation, including the VME. 
The values, g pnl = 0.103 and gr WW7 = 0.313 are adopted. 

In the 7r + production, contributions from the VME are found to be negligible. One reason 
is the relatively smaller couplings of g piT1 and g P NN compared to the couplings for u exchange. 
However, the main factor leading to small VME contributions in the tc + production is a large 
cancellation occurring among terms proportional to Qa and g t in Eq. pBj). Since different 
contributions from the VME to different reactions depend on the quark model prediction for 
gA and gt, the VME might introduce more model-dependent ingredients in the calculations. 



IV. CONCLUSIONS 

We have studied pion photoproduction in four charge channels within the quark model 
incorporating an effective Lagrangian. Up to E 1 = 500 MeV, the cross sections and single 
polarization asymmetries can be accounted for with one adjustable parameter for the Delta 
excitation strength. We find that if a stronger coupling for the Delta transition is employed, 
all the observables in the first resonance region can be reproduced. In another words, the 
NRCQM might have provided a reasonable form factor for the Delta resonance, but with 
weaker coupling. 

As the first systematic microscopic study of pion photoproduction, this result suggests 
that a direct calculation of the tree level diagrams based on the quark model with a chiral 
effective Lagrangian may contain the main ingredients required in an elementary approach. 
In particular, we show that the quark model form factors play a key role in reproducing 
the data over a wide energy region. Such a form factor can be only self-consistently and 
completely considered in a direct calculation of quark level diagrams. This result highlights 
the relation between the background Born terms and resonance excitations. Extensions to 
higher energies would help to clarify the relation between quark model results and those 
found via phenomenology. Nevertheless, a parallel investigation of electroproduction would 
be useful for a better understanding of the Delta resonance based on this model. We shall 
report this elsewhere. 

Restricted to the low energies at E 1 < 500 MeV, we see that t-channel vector-meson 
exchange is negligible. This leaves the duality hypothesis far from conclusive. We expect 
that more sophisticated calculations at high energies could be helpful in disentangling this 
question as well. 
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FIG. 1. Tree level diagrams calculated in this model. 
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FIG. 2. Real and imaginary parts of the multipole M l+ derived from the Delta amplitude in 
7T° channel. Model results with and without the quark model form factor (FF) are compared. Data 



are from SAID analysis [27]. 
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FIG. 3. Differential cross sections for 7p — * ir + n. The solid curves denote full calculations with 
Ca = 1.7474; Dotted curves, results in the SU(6)®0(3) symmetry limit; Dashed curves, results 
without n > 1 it-channel resonance contributions; Dot-dashed curves, calculations with resonance 



real parts eliminated. Data are from Ref. [28-3 
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FIG. 4. Polarized beam asymmetry E for ^yp — * ir + n. The solid curves denote full calculations, 
while the dashed present results without n > 1 it-channel resonance contributions. The dotted 
curve at = 260 MeV denotes the effects by eliminating the nucleon pole terms. Data are from 
Ref. pp]. 



23 



h 1 

0.8 
0.6 
0.4 
0.2 




E =220 MeV 




50 100 150 




50 


100 150 


~ / i i i i 1 i i i i 1 i i 


45 MeV 

1 I I \ 


50 


100 150 


\/$ 

l/l 1 ! ! 1 ! 1 1 


w 

iii 


91 MeV 


50 


100 150 


i i i i i i i i i 


E T =700 MeV 
«^*> 

I I I T T T 1 I I 



50 100 150 
6 (deg) 



50 100 150 
6 (deg) 



FIG. 5. Polarized target asymmetry T for — > ir + n. Data are from Ref. . 41,49— 51 1 . 
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FIG. 6. Recoil polarization asymmetry P for —* 7r + n. Data are from Ref. [ 41 , 52 1 
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FIG. 8. Differential cross sections for ~/p — ► 7r°p. Data are from Ref. HHH 
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FIG. 10. Polarized target asymmetry T for ^yp — > 7r°p. Data are from Ref. [70,75-77,74] 
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FIG. 11. Recoil polarization asymmetry P for — > TT°p. Data are from Ref. [78,70,72,74]] 
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FIG. 12. Differential cross sections for -yn — ► tt n. Data are from Ref. |8C] 
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FIG. 13. Cross sections for ^yp — > 7r + n at 300 MeV. Plotted are the Born terms with (solid 
curve) and without (dotted) form factor, and for the Born terms plus the Delta transition with 
(dashed) and without (dot-dashed) form factors, respectively. 
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FIG. 14. Energy evolution of the Born terms plus Delta helicity amplitudes compared with the 
SAID analyses. Columns from left to right are amplitudes ^1,2,3,4- The solid and dashed curves 
denote the real and imaginary part calculated by this model, while the dotted and dot-dashed 
denote those by SAID analyses. 
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FIG. 15. Observables for jp — > 7r + n at three energies. Solid curves denote results for s+u; 
dashed for s+u+t; and dotted for s+t Data are from Ref. . 
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